G\ 1 Nuclear constraints on the inner edge of neutron star crusts 
O ' 



o 

a Department of Physics, Shanghai Jiao Tong University, Shanghai 200240, China 

O 

m 



=5 



c3 



Lie- Wen Chen a *, Bao-An Li b t, Hong-Ru Ma a , and Jun Xu a| 



b Department of Physics and Astronomy, Texas A&M University- Commerce, Commerce, 
Texas 75429-3011, USA 



1. INTRODUCTION 



We show that the widely used parabolic approximation to the Equation of State (EOS) 
of asymmetric nuclear matter leads systematically to significantly higher core-crust tran- 
sition densities and pressures. Using an EOS for neutron-rich nuclear matter constrained 
by the isospin diffusion data from heavy-ion reactions in the same sub-saturation density 
range as the neutron star crust, the density and pressure at the inner edge separating 
the liquid core from the solid crust of neutron stars are determined to be 0.040 fm -3 
< Pt < 0.065 fm~ 3 and 0.01 MeV/fm 3 < P t < 0.26 MeV/fm 3 , respectively. Implications 
■ of these constraints on the Vela pulsar are discussed. 

in 

On 

The inner crust in neutron stars spans the region from the neutron drip-out point to 
the inner edge separating the solid crust from the homogeneous liquid core and plays an 
important role in understanding many astrophysical observations [ [U El E] • While the 
neutron drip-out density p out is relatively well determined to be about 4 x 10 11 g/cm 3 [ 
H], the transition density p t at the inner edge is still largely uncertain mainly because of 
our very limited knowledge on the EOS of neutron-rich nucleonic matter, especially the 
density dependence of the nuclear symmetry energy E sym (p) [[2]. 

Recently, significant progress has been made in constraining the EOS of neutron-rich 
nuclear matter using terrestrial laboratory experiments (See, e.g., Ref. [|S] for the most 
recent review) . In particular, the analysis of isospin-diffusion data [ El [71 E] in heavy- ion 
collisions has constrained tightly the E sym (p) in exactly the same sub-saturation density 
region around the expected inner edge of neutron star crust. In the present talk, we report 
our recent work on locating the inner edge of the neutron star crust using terrestrial nuclear 
laboratory data within both thermodynamical and dynamical methods [E]. 
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Figure 1. (Color online) The p% verus L from 
the dynamical and thermodynamical methods 
with and without the PA using the MDI in- 
teraction. The triangles are obtained by Ku- 
bis [[11] and the star with error bar represents 
L = 86 ± 25 MeV. Taken from Ref. [[9] 



Figure 2. (Color online) The Pt versus L 
and pt by using the dynamical method with 
and without the PA using the MDI interaction. 
The star with error bar in the left panel rep- 
resents L = 86 ± 25 MeV. Taken from Ref. [ 
[9] 



2. RESULTS 

Shown in Fig. [1] is the pt as a function of the slope parameter of the symmetry energy 
L = 3po dEsy d ™ ip) \p= PO with the MDI interaction [[10]. For comparisons, we have included 
results using both the dynamical and thermodynamical methods with the full EOS and its 
parabolic approximation (PA), i.e., E(p,5) = E(p,5 = 0) + E sym (p)5 2 + 0(5 A ) using the 
same MDI interaction. The thermodynamical method is the long wave length limit of the 
dynamical one when the Coloumb and surface interactions are neglected [19]. With the full 
MDI EOS, it is clearly seen that the p t decreases almost linearly with increasing L within 
both methods. It is interesting to see that the two methods give very similar results (the 
difference is actually less than 0.01 fm~ 3 ). Surprisingly, the PA drastically changes the 
results, especially for stiffer symmetry energies (larger L values). Also included in Fig. [T] 
are the predictions by Kubis using the PA of the MDI EOS in the thermodynamical 
approach [IT]. The large error introduced by the PA is understandable since the /in- 
stable npe matter is usually highly neutron-rich, thus the contribution from the higher 
order terms in 5 is appreciable. This is especially true for the stiffer symmetry energy 
which generally leads to a more neutron-rich npe matter at subsaturation densities. In 
addition, simply because of the energy curvatures involved in the stability conditions, 
the contributions from higher order terms in the EOS are multiplied by a larger factor 
than the quadratic term. Our results indicate that one may introduce a huge error by 
assuming a priori that the EOS is parabolic with respect to the isospin asymmetry for 
a given interaction in calculating the pt. We thus apply the experimentally constrained 
L to the p t — L correlation obtained using the full EOS in constraining the p t in the 
dynamical method as it is more complete and realistic. As shown in Fig. [TJ the constrained 
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L = 86 ±25 MeV obtained from the transport model analysis of the isospin diffusion data [ 
El El E] then limits the transition density to 0.040 fm -3 < p t < 0.065 fm -3 . 

The pressure at the inner edge, P t , is also an important quantity which might be 
measurable indirectly from observations of pulsar glitches [|2j [3]. Shown in Fig. [2] is the 
P t versus L and p t . Again, it is seen that the PA leads to huge errors for larger (smaller) 
L (p t ) values. For the full MDI EOS, the P t decreases (increases) with the increasing L 
(pt) while it displays a complex relation with L or p t for the PA. The complex behaviors 
are due to the fact that the p t does not vary monotonically with L for the PA as shown 
in Fig. [TJ From the constrained L values, the Pt is limited between 0.01 MeV/fm 3 and 

0. 26 MeV/fm 3 . 

The constrained values of pt and Pt have important implications on many properties 
of neutron stars [HI [21 [3]. As it was shown in Ref. [[2], the crustal fraction of the total 
moment of inertia of a neutron star, i.e., A I / 1 depends sensitively on the P t and p t at 
subsaturation densities, but there is no explicit dependence upon the higher- density EOS. 
So far, the only known limit of AI/I > 0.014 was extracted from studying the glitches of 
the Vela pulsar [ 3J. This together with the upper bounds on the P t and p t (pt = 0.065 fm -3 
and P t = 0.26 MeV/fm 3 ) sets approximately a minimum radius of R > 4.7 + 4.OM/M 
km for the Vela pulsar. The radius of the Vela pulsar is predicted to exceed 10.5 km 
should it have a mass of 1.4M . A more restrictive constraint will be obtained from 
the lower bounds of p t = 0.040 fm~ 3 (P t = 0.01 MeV/fm 3 ) and it can be approximately 
parameterized by R = 5.5 + 14.5M/M Q km. It is thus seen that the errors in both the 
transition density and pressure are still large. Thus, the uncertainties for the mass-radius 
relation of the Vela pulsar are still large. A conservative constraint of R > 4.7 + 4.OM/M 
km using the upper bounds on the P t and p t was then obtained [ [H] • We notice that a 
constraint of R > 3.6 + 3.9M/M km for this pulsar was previously derived in Ref. [[3] by 
using p t = 0.075 fm -3 and Pt = 0.65 MeV/fm 3 . However, the constraint obtained using 
data from both the terrestrial laboratory experiments and astrophysical observations is 
significantly more stringent. 
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